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Daniele Vilone, 1 Alessandro Vespignani, 2 and Claudio Castellano 1 ' 3 

1 Dipartimento di Fisica, Universita di Roma "La Sapienza" , P.le A. Moro 2, 1-00185 Roma, Italy 
2 The Abdus Salam International Centre for Theoretical Physics (ICTP), P.O. Box 586, 34100 Trieste, Italy 
C<\ I S INFM, Unitd di Roma 1, P.le A. Moro 2, 1-00185 Roma, Italy 

O (Dated: February 1, 2008) 

O ■ 

^s^j . We study the one-dimensional behavior of a cellular automaton aimed at the description of the 

formation and evolution of cultural domains. The model exhibits a non-equilibrium transition 
' between a phase with all the system sharing the same culture and a disordered phase of coexisting 

("S , regions with different cultural features. Depending on the initial distribution of the disorder the 

transition occurs at different values of the model parameters. This phenomenology is qualitatively 
captured by a mean-field approach, which maps the dynamics into a multi-species reaction-diffusion 
problem. 
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PACS numbers: 87.23.Ge, 05.50.+q, 89.75.Fb 



I. INTRODUCTION 



In recent years cellular automata and lattice models have been introduced in many fields in order to investigate 
the emergence of nontrivial collective phenomena out of simple microscopic interactions between a large number of 
elementary objects jj], |2|, |3|. In the context of social phenomena, numerical simulations of simple automata have been 
used to study the formation and spreading of cultural domains I | a |, Culture is intended here as the set of 
behaviors, beliefs, technical standards or values that individuals may mutually exchange. The cultural features are 
modeled as a set of F variables, which can assume integer values between 1 and Q, specifying the various traits that 
each feature can assume [Q. In general, the environment where exchanges take place is a two-dimensional square 
grid, where each site represents an individual or a group of individuals, whose culture is defined by the specific values 
of the F variables. Starting from an initial random configuration, interactions take place between nearest neighbors, 
leading to the formation and evolution of regions of homogeneous culture. It is obvious that the assessment of how 
realistic the model is or to which extent its results can be useful in the interpretation of empirical data in the social 
science context is a task left to sociologists. From the point of view of statistical physics this model (Axelrod model) 
is appealing, because it has a non-trivial out of equilibrium dynamics, similar but not equivalent to other well studied 

■ models §,§. 

In previous works ^, ||, |[ |t| , it has been investigated how the model behavior changes with respect to the values of 

■ the parameters F and Q. In particular it has been shown by large-scale simulations that, on a two-dimensional lattice, 
the model exhibits a non-equilibrium transition between a phase with all the system sharing a common culture and 
a disordered phase of coexisting regions with different cultures. This is a very interesting example of ordering phase 
transition in the context of far from equilibrium statistical physics. In this paper we perform a thorough analysis of 
the model defined in Ref. on a one-dimensional lattice. This geometry acquires a special role in the description 
of social systems, since it can be used as starting point for complex topologies such as small- world networks |L0, pd| . 
Noticeably, in dimension d = 1 it is possible to provide an analytical treatment of the model by mapping the evolution 

i-^J ■ dynamics into a multi-species reaction-diffusion system, that we study within a mean-field approximation. We recover 
also in one dimension the existence of an ordering phase transition, occurring at a critical value of the parameter Q. 
When the initial value of the cultural variables on each site is chosen randomly with uniform distribution between 1 
and Q, we find that the system undergoes a transition between a fully ordered state for Q < Q c and a disordered one 
J> , for Q > Q c . When the distribution of initial variables is instead Poissonian of parameter q, the system exhibits a phase 
transition between order and disorder, but only for F > 2. For F — 2 any finite value of q leads, in the thermodynamic 
limit, to an asymptotic disordered state. In order to test the analytical predictions we perform numerical simulations 

■ of the original dynamics, finding a good agreement with the theoretical results. 
The paper is organized as follows. In Section [n] the model is introduced, the dynamics qualitatively described and 



the results of the two-dimensional case briefly summarized. In Sec. Ill we present the mapping to the reaction-diffusion 
system and the mean-field approach to it. In Sec. |y| we report the results of numerical simulations of the model in 
one dimension, both for the case of uniform and Poisson initial distributions, and compare them with the analytical 
predictions. Finally in Sec. [v| we draw some conclusions. 
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II. THE MODEL 

We consider a one-dimensional lattice of linear size L. On each site i there are F integer variables <Jij which define 
the cultural features of the individuals living on that site. In the original model 0, each feature / = 1, . . . , F on each 
site i is initially drawn randomly from a uniform distribution on the integers between 1 and Q. The parameter Q is 
a measure of the initial cultural variability (i.e. disorder) in the system. In Ref. |(| instead, the initial values were 
chosen according to a discrete Poisson distribution of parameter q Prob(<7ij = k) = q k e~ q /k\, so that the positive 
real q was the average of the values extracted. Results were found to be qualitatively similar independently from 
the initial distribution, and only the Poissonian case was studied in detail. Here we will consider also the case with 
uniform distribution and we will discuss analogies and difference between the two cases. 

The dynamics is defined as follows. At each time step, a pair of nearest neighbor sites i and j is randomly chosen. 
A feature / is selected and if Oij ^ <Zj,/ nothing happens. If instead <Jij — <jjj then an additional feature /' is 
randomly chosen among those taking different values across the bond, Uiji ^ ®j.f- Such a feature is then set equal 
in the two sites: o% H j> — > cr^ fi = Vj.f- Time is measured as the total number of attempts divided by the number 
of sites L. Axelrod model can be seen as F coupled voter models The evolution of the system is basically a 
coarsening dynamics of ordered regions. The initial state is totally disordered since the choice of variables in each 
site is uncorrelated to any other one. As the dynamics starts, the interaction between sites tends to make neighbors 
similar, introducing correlations in the system. Ordered regions start to form with all sites sharing exactly the same 
value of all variables. Clearly there are many different possible types of ordered regions. As time proceeds boundaries 
between regions move, and the average domain size increases. On a finite lattice the ordering process continues until 
the system reaches an absorbing state such that no further dynamics is possible. Absorbing states can be of two types. 
A perfectly ordered absorbing state is made by a single region covering the whole system: all sites have the same set of 
variables cr%j- Otherwise an absorbing state can also be made by different ordered regions such that adjacent regions 
have all variables different. In that case all sites are either totally equal or totally different from their neighbors and 
no interaction may take place. 

The nature of the final absorbing state depends on the full evolution of the system, i. e. it is the result of the 
competition between the disorder of the initial configuration and the drive toward order due to local interactions. 
Intuitively one expects that for small initial disorder (small Q or q) local interactions will prevail and the final state 
will be perfectly uniform, while in the opposite limit the ordering process will not last long enough to fully overcome 
the disorder of the initial condition. This picture turns out to be correct in two dimensions filj. In that case extensive 
numerical simulations have shown that for small values of q, i. e. small initial disorder, the model converges toward 
an ordered state where a single culture occupies a macroscopically large fraction of the system. On the contrary 
for large values of q many different cultures coexist in the final absorbing state. In the thermodynamic limit of an 
infinitely large system one sees a well defined non-equilibrium phase transition for a precise value q c , separating an 
ordered phase (for q < q c ) from a disordered one (for q > q c ). Interestingly, the nature of the transition is different 
depending on the number F of variables used to define the culture of a site. If F = 2 the transition is continuous, 
with the order parameter (defined below) vanishing for q — > q c . For F > 2 instead the order parameter exhibits a 
jump at the transition. The different behavior is reflected also in the distribution of region sizes at the transition. 
For F = 2 it is a power law with an exponent r < 2, while for F > 2 the exponent is universal and larger than 2, so 
that the transition does not involve the divergence of a correlation length Q . It is important to stress that the kind 
of nonequilibrium phase transition discussed here is not an absorbing phase transition jl| . Here the system always 
reaches a frozen absorbing state. The transition occurs in the properties of such absorbing state and hence is static 
in nature, although it is clearly determined by the dynamical evolution of the problem. 

In this paper we shall investigate whether the same kind of transition occurs also in a one-dimensional geometry. 
To monitor the evolution of the model we mainly focus on the density (i. e. the number divided by the system size 
L) of active bonds n>A(t) and of fixed bonds np(t). We define as active a bond between two adjacent sites such that 
they have more than and less than F equal features. We define as fixed a bond between two adjacent sites such 
that they have no equal variables. Clearly, active bonds are where the dynamics take place and form the boundary 
between regions that may interact, leading to further ordering. Boundaries formed by fixed bonds are instead frozen. 
All bonds connecting sites belonging to the same region are neither active nor fixed. We denote their density as Uq. 
Clearly, for all times no + ha + np = 1. Both the density ua of active bonds and the density np of fixed bonds 
measure correlations in the system: 1/np and 1/ha are correlation lengths growing in time. Since for any finite lattice 
an absorbing state is always reached by the dynamics, n,A(t) goes asymptotically to zero in all cases. The value of 
np{t = oo) measures the order in the system. If np(oo) keeps finite as L diverges, the system remains disordered in 
the thermodynamic limit, while np{oo) — > indicates full order. 

Another useful quantity that will be considered below, is (S max ) / L, where {S ma x) is the average over different 
realizations of the size of the largest region existing in the final frozen state ||. {S max )/L is the fraction of the whole 
lattice occupied by the largest region and is an alternative order parameter, being zero for a disordered state and 
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FIG. 1: Example of the evolution of three neighboring sites for F — 4. Site (B) and its neighbor (A) are selected and their 
feature / = 1 is compared. Since <ta,i = crs.i = 5, the two sites interact and feature /' = 4 of site (B) is set equal to ctaa = 4. 
Bond (A)-(B), which was of type 3 becomes of type 2. Bond (B)-(C), initially of type 2 becomes of type 3. 

finite when a macroscopically large region exists in the system. 

In the next Section we describe an analytical treatment of the problem based on the mapping of the dynamics into 
a multi-species reaction-diffusion system, which is then analyzed in terms of mean-field equations. This approach 
turns out to capture remarkably well the qualitative behavior of the system. In some cases it gives predictions that 
are quantitatively accurate. 

III. THEORETICAL APPROACH 
A. Mapping to a reaction-diffusion system 

The Axelrod model considers the sites as elementary objects, which interact according to the rules described above. 
This dynamics can also be thought in terms of bonds between sites whose dynamics is induced by the updating of 
sites. More precisely, let us define as bond of type k a bond connecting two sites with k different features. Clearly k 
can be any number between and F. Bonds of type connect sites with all variables equal, i. e. belonging to the 
same region. Bonds of type F instead connect sites completely different and are therefore the fixed bonds introduced 
above. Bonds of type 1, . . . , F — 1 are what we called active bonds in the previous section. 

To understand how the update of sites induces a dynamics for bonds, let us consider three neighbor sites (A), (B), 
(C). Site (B) is updated depending on the configuration of its neighbor (A). If the bond between sites (A) and (B) 
is of type or F nothing occurs. Otherwise, if the bond is of type n one variable on site (B) is set equal to the 
corresponding variable on site (A). Hence the bond (A)-(B) becomes of type n — 1. The update of site (B) modifies in 
general also the state of bond (B)-(C), which initially was of type m. Notice that for this bond the number of features 
which arc different may increase, decrease or even remain unchanged, so that after the update the bond may be of 
type m — 1, morm+1. In summary, if we denote a bond of type k as Ak, the update of a site is then equivalent to 
a two-particle reaction 

A n + A m > A n —i + A m +j C? = o,±i). (f) 

A specific example is presented in Fig. [|. Bonds of type can be considered as a sort of vacuum state where other 
types of bonds move. Reactions of the type A\ + Aq — > Ao + A\ are then diffusion processes for bonds of type 1. 
These are the only diffusion processes allowed. 

In order to fully specify the dynamics we must provide the transition rates for the reactions (|f|). Let us consider 
the three sites of Fig. ||. One of the features of site (B) is selected and compared with the corresponding one 
of site (A). In order for a reaction to occur the two variables must be the same. This occurs with a probability 
Rn^ = (1 — n/F)(l — 8 n ,o) where we have taken into account that nothing happens if the bond is of type 0. 

Hence with probability Rn \ another variable <JbJ' is set equal to <JAJ'- What is the effect of this change on the 
bond (B)-(C)? If crc,f wa s equal to <7bj', then the update increases by one the number of features different in the 
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bond (B)-(C). The probability that crcj' — °~B.f' is the probability that we randomly select one of the F — m equal 
variables out of the F total variables 
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With probability 1 — R m _ x instead, the variable selected will be different in site (B) and (C), o"c,/' G Bj'- In this 
case, since <JbJ' is set equal to caj' the final state of bond (B)-(C) depends on whether <jaj' and o~c,p are equal or 
not. Let us denote as A the probability that <taj' — °~C,f- Then the probability that bond (B)-(C) goes from type 
m to type m — 1 is 



d(2) 
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The probability R rn o i s obtained via the normalization condition 
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In summary the reaction A n + A m A n —\ + A m+ j occurs with rate 



(5) 



The reactions (|l|) and the rates (||) fully specify the dynamics of the model. Provided one is able to write down the 
appropriate form of A the mapping from the original Axelrod model is exact and Eqs. (|l|) and (^|) are an alternative 
formulation of the same dynamics. 



B. Mean- field approximation 

As discussed above, A is the probability that two independent random variables (A and C), different from a third 
one (B), are equal. Here we will determine it by introducing a mean-field approximation, i. e. assuming spatial 
uniformity. 

Let us call k the value of variable (B) and ps(fc) its distribution. The A can be written as ^2^PB{k)\k where A& is 
the probability that A and C assume the same value, different from k. Clearly 

Afe = PA(i)Pc(j)Si,j (6) 

where pA(i) — PA{i)/ Si^fe v(i) is the probability distribution of the variable A constrained to be different from k. 
We assume that the distributions of all variables A, B and C remain equal to the initial ones at all times. For uniform 
initial distribution p(i) — 1/Q, we recover the expected value 

(7) 



Q-i 

For variables distributed according to a Poissonian, the plot of A as a function of q, is presented in Fig. |^. The two 
behaviors are qualitatively similar. However, one has to remark that in the Poissonian case A = 1 for q = 0, in the 
trivial situation when all variables are the same. For uniform initial condition instead A = 1 for Q = 2, i. e. when 
the variables may assume two values. 

From the knowledge of the possible reactions and their associated transition rates it is straightforward to write 
down the mean-field equations for the densities rik of bonds of type k. 

We obtain the set of equations {k = 0, . . . , F) 

n k = -R^rik + R^nk+i (8) 

F-l 

+ £ R$n N [-n k (l - + n k -iR<V ltl + n^R™^ 
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FIG. 2: Plot of A as a function of Q (uniform case) or q (Poisson case) 



where clearly n_i = = np + i. More explicitly 

n k = - ( 1 - 4 J (1 - 4,oW + f 1 - — ^ ) ^fc+i (9) 
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It is easily verified that X)fc=o ™ k = 0. Eqs. (^|) are mean-field equations, since they neglect spatial fluctuations and 
noise in the system. They are complemented by the initial conditions, which are 

n k (t = 0) = pF- k (l-P eq ) k (^ (10) 

where P eq is the probability that a feature takes the same value in two adjacent sites. Clearly P eq — 1/Q in the 
uniform case, while 

h=0 y ' 

for a Poisson initial distribution. 

The mean-field system of equations of motion (||) consists of F coupled nonlinear equations. A detailed study of the 
temporal evolution is necessarily performed by integrating them numerically. However some important conclusions 
on the asymptotic state of the system can be drawn analytically by considering the fixed points of Eqs. @. A first 
crucial observation is that, for any F, a state with all = (k = 1, . . . , F — 1) is a fixed point of the equations. 
Therefore the system always admits a family of frozen asymptotic states, parametrized by the density of fixed bonds 
n p (no = 1 — riF because of normalization). To understand which of these fixed points is selected by the dynamics 
we perform a stability analysis, by linearizing the dynamical system around the generic fixed point = 0. 

The stability matrix has F eigenvalues. One eigenvalue is always zero, corresponding to the direction parametrized 
by np and all made of fixed points. F — 2 other eigenvalues are real negative numbers or complex ones with negative 
real part, corresponding to stable directions. The last one, that we denote as is always real and it crucially depends 
on the value of A Its value as a function of A, computed analytically for F = 2 and F = 3, is plotted in Fig. || for 
two different values of np. Clearly [i depends from the value of rip, but for any np > it is negative for A < \/{F — 1) 
and positive for A > l/(F — 1). For np = (homogeneous absorbing state) /i = for any value of A, so that the 
stability of this fixed point cannot be assessed directly within a linear analysis. However its stability can be inferred 
indirectly from the stability of other fixed points: if all other fixed points are unstable the fixed point with np = 
must be stable; if all others are stable, the dynamical system will not go to np = 0. This allows us to conclude that 
for values of the parameter Q or q of the model such that A < 1/(^—1) the equations will converge exponentially 
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FIG. 3: Plot of the eigenvalue jj, versus A for F = 2 (thin lines) and F = 3 (thick lines), and uf = 0.5 (solid) and uf — 0.8 
(dashed). 
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FIG. 4: Mean-field phase diagram of the model. The two critical lines are for uniform initial distribution (circles) and for 
Poisson initial distribution (diamonds). 



fast to a finite value np > implying that the system remains disordered. If instead A > 1/(F — 1) the only possible 
attractive fixed point is uf = 0, which will attract the solution more slowly than exponentially. In this case the 
system exhibits a power-law convergence to an ordered state. Whether this transition actually takes place and for 
which value of the parameter of the initial distribution, it depends on the precise value of A, that is plotted in Fig. |^. 
For uniform initial distribution, since A = 1/(Q — 1) it is immediate to see that a phase transition occurs for 



= F 



(12) 



so that for Q < Q c the system is ordered, while the system remains asymptotically disordered for Q > Q c . 

For Poisson initial distribution and F > 2, the scenario is the same, with an ordered phase for small q separated 
by a sharp transition from an ordered phase for large q. For F = 2 all finite values of q imply a positive eigenvalue 
/i and hence a disordered state. Only for q = order is recovered, but of a very trivial type, since q = means that 
all variables are the same already in the initial state. Notice that, since A decays as q" 1 ^ 2 for large q, the transition 
point diverges for large F as q c ~ F 2 , differently from the uniform case. 

For easier insight into the behavior of the model we plot in Fig. ^ the mean-field phase diagram. 
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C. Uniform initial distribution with F — 2 



Let us now give a closer look at the case of uniform initial distribution with F = 2. In this case the set of mean-field 
equations ^ is reduced to two coupled equations 



Til 
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(13) 
(14) 



where no has been eliminated using normalization. In the case Q — 2 it is easy to solve (|T4| ) explicitly, since the 
equation for ni does not depend onnj?. One obtains 



m(t) = 

n F (t) 



1 



(2 + i)V2 2(2 + «)' 



(15) 
(16) 



Hence we find, as expected, that both the density of active bonds and of fixed bonds vanish asymptotically as 
power-laws. 

For Q > 2 the two equations (|l4|) become fully coupled. The asymptotic value of np is found to be 



np(t = oo) 



Q-2 
Q 



(17) 



for all values of Q, which will be shown below to agree very well with numerical results. 



IV. NUMERICAL RESULTS 



In order to validate the results of the theoretical approach we compare them with extensive numerical simulations 
of the model. For each value of the parameters, we average over 10000 or more realizations. For both types of initial 
distributions, all values of F and of q or Q, we find that the system reaches for long times an absorbing state with no 
active bonds (ua — 0). Extrapolation to the limit L — > oo indicates that also in the thermodynamic limit the system 
freezes asymptotically. This is the first prediction of the mean-field approach that is verified numerically. We now 
turn to a more detailed analysis of the asymptotic state as parameters are changed. 



A. Uniform initial distribution 



Let us consider first the case with F = 2 variables on each site, which in two dimensions gives a continuous 
transition. The mean-field approach predicts the existence of a transition between an ordered and a disordered phase 
for Q c = 2. Such a prediction is confirmed by the behavior of the order parameter (S max ) / L. As Fig. || illustrates, 
(Smax)/L keeps for Q = 2 a finite value close to 0.71 in the thermodynamic limit, while it goes to zero for Q > 2. 
Hence one macroscopic region occupies a finite fraction of the system for Q = 2, while for larger values of Q, no 
macroscopic domain is formed. 

We remark that in the ordered case the asymptotic value of (S ma x)/L is not 1. This should not be taken as evidence 
that the system reaches a partially ordered state with the largest region spanning only a finite part of the whole lattice. 
In fact the situation is more complex. (S max ) is the average over a large number of realizations. In many of them the 
system becomes fully ordered and Smax/L = L However in a finite fraction of the realizations the system remains 
partially disordered so that the order parameter has a value smaller than 1. In practice the distribution of the values 
of S m ax/L for L — > oo is a ^-function for S max /L = 1 plus a finite nontrivial contribution for S max / L < 1 This leads 
to (S max )/L ~ 0.71. 

For Q > 2 the system remains disordered and a finite fraction np of the bonds is fixed. In the previous Section, 
Eq. ( |l7j ) is the mean- field result for np in the asymptotic state as a function of Q. Such a formula is compared 
with numerical results in Fig. ^. The agreement is excellent, indicating that the mean- field approach is quantitatively 
correct in this case. 

Let us now analyze the whole dynamical evolution. In the main part of Fig. [?] we plot the temporal behavior of the 
density of active bonds in the case Q = 2, for different system sizes. All curves collapse on the same one when 
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FIG. 5: The order parameter (S m ax)/L as a function of the system size L for F = 2. From top to bottom the curves are for 
Q = 2, Q = 3, Q = 4 and Q = 5. 
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FIG. 6: Comparison of the numerical results with the mean-field formula (Q — 2)/Q for the final density of fixed bonds in the 
case F — 2. 

time is rescaled by L 2 , i. e. n A obeys the scaling form 

n A (t,L)=t-^ 2 g A (t/L 2 ). (18) 

The scaling function g A (x) goes to a constant as x <C 1 and goes to zero exponentially for i» 1. In the inset of the 
same figure it is shown that also the density of fixed bonds obeys a scaling form 

n F (t,L)=t- 1 / 3 g F (t/L 2 ) (19) 

where gF(x) — const for x <C 1 and gF(x) ~ x 1 ^ 3 for x 3> 1. 

The physics behind Fig. is clear. In the initial state a finite number of active bonds is present. Notice that 
for F = 2 and Q = 2 a single type of active bond exists, connecting sites with one equal and one different feature. 
Interaction between sites implies diffusion of active bonds and annihilation whenever two active bonds meet. Fixed 
bonds are affected by active bonds but the converse is not true: As far as active bonds are concerned, the system is 
equivalent to the reaction-diffusion system A + A — > 0. The decay as t -1 / 2 of the density of active bonds is easily 
related to the combination of diffusive wandering and mutual annihilation |13j . The inverse of n A is a correlation 
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FIG. 7: Main: Scaling plot of the density of active bonds ua versus time for F = 2, Q = 2 and system size (left to right) 
L = 10, 20, 50, 100, 200, 400, 1000, 2000, 5000. Inset: Scaling plot of the density of fixed bonds n F versus time for the same set 
of parameters. 

length that grows in time. In an infinite system such length grows indefinitely. In a finite lattice, growth stops when 
X/ua reaches the system size L, and this yields a cutoff time scaling as L 2 . The evolution of fixed bonds is enslaved 
by the dynamics active bonds but it is slower, l/np ~ t 1 / 3 is a correlation length non equivalent to 1/ua ~ t x l 2 . It 
grows with time but remains finite when the cutoff time is reached, so that the in the frozen state np reaches a finite 
value decreasing as L -2 / 3 , which vanishes in the thermodynamic limit. 

A comparison of the exponents found numerically with the mean-field values [Eq. indicates that the mean-field 
approximation correctly reproduces that uf decays more slowly than tia, but it does not predict the correct exponents. 
This is not a surprising result: we have computed the mean-field rates assuming the system to be homogeneous, while 
instead ordered regions grow indefinitely. The problem is here the same of the mean-field treatment for the reaction- 
diffusion system A + A — > |l|. For higher values of Q the disagreement between the mean-field and numerical 
results is similar. It is somewhat striking that despite the dynamics is not captured quantitatively by the mean-field 
approximation, the asymptotic value for np is remarkably accurate (Fig. ^|). 

Let us now discuss the case F > 2 corresponding to a discontinuous transition in d = 2. For illustrating the behavior 
common to all values F > 2 let us focus on the case F — 10. As in the two-dimensional case ||, the system exhibits 
a discontinuous transition between a fully ordered phase for small Q and a disordered phase for Q > Q c . This can be 
inferred from Fig. ||, where 1 — {S max )/L is plotted versus Q. By observing that the convergence of {S max )/L to 1 is 
exponential for Q < 6, while curves are bent upward for Q > 7 one is lead to conclude that the transition occurs for 
Qc * 7. 



B. Poissonian initial distribution 

Let us consider now a system such that the initial values of the variables are chosen according to a discrete Poisson 
distribution of parameter q. q plays here a role similar to the parameter Q for the uniform initial distribution, since it 
fixes the average value of the variables in the initial state. However, no matter how small is q, arbitrarily large values 
can be extracted with a Poisson distribution. 

Let us consider first the case F — 2 (Fig. ^|). The behavior differs qualitatively from the uniform case. For any 
value q > it is clear that the order parameter goes asymptotically to zero, even if for small q, (S max )/L remains 
close to 1 for relatively large L. We conclude that there is strong evidence that for F = 2 an initial distribution of 
Poisson type destroys long range order for any value of q > 0, as predicted by the mean-field approach. 

We finally turn to the case F = 10, as an example of what happens in general for F > 2. The value of 1 — (S max ) /L 
is plotted in Fig. [l(]. For small q it is clear that the system becomes ordered as L grows while for large q, the curves 
tend to go toward 1, indicating the existence of a discontinuous phase transition for intermediate q. A rough estimate 
of the transition point is between q — 4 and q = 5. 

These numerical data corroborate the picture emerging from the mean-field approach. For F = 2 we actually find 
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FIG. 8: Plot of 1 — {S ma x)/L as a function of the system size L for F — 10. From bottom to top the curves are for Q — 3, 
Q = 5, Q = 6, Q = 7, Q = 8, Q = 9, Q = 10 and Q = 15. 
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FIG. 9: Plot of 1 — {Smax)/L as a function of the system size L for F — 2 for Poisson initial distribution. From bottom to top 
the curves are for g = 0.02, q — 0.05, g = 0.1, q = 0.25, q = 0.5, q = 1, g = 2 and g = 10. 

no evidence of long range order in the system, no matter how small q is taken. For F = 10 we observe instead the 
pattern already found with uniform initial distribution: a transition between a disordered phase for large q and an 
ordered one for smaller values of the parameter q. Also the location of the transition seems to be captured with good 
accuracy. 

V. CONCLUSIONS 

In this paper we have analyzed in detail the behavior in one dimension of the Axelrod model, a cellular automaton 
introduced in Ref. Q for the description of social influence. Interestingly, the one-dimensional version of the model 
displays a phenomenology as rich as the two-dimensional case. More precisely, we find by numerical simulations 
a non-equilibrium phase transition between a disordered and an ordered phase. This phenomenology is recovered 
accurately by a simple mean-field treatment of the problem, which exploits a mapping onto a reaction-diffusion 
problem involving F different species. The mean-field approach is successful in several aspects. First, it correctly 
predicts that the system always reaches an asymptotically absorbing state. Secondly, it captures the existence of a 




FIG. 10: Plot of 1 — (Smax) I L as a function of the system size L for F — 10 for Poisson initial distribution. From bottom to 
top the curves are for q = 1, q = 2, q = 3, q = 4, q = 5, q = 6, q = 8, and q = 10. 



transition in the final absorbing state or the absence of it (for Poissonian initial conditions and F = 2). Finally, it 
provides reasonably accurate estimates of the critical value of the control parameter Q or q. Mean-field fails instead 
in the prediction of the exponents of the temporal behaviors close to the transition. This is a typical scenario for 
mean-field treatments: transition points can be computed with good accuracy while exponents are not captured 
correctly. The theoretical approach presented in this paper is the simplest suitable for this kind of problem. For more 
precise values of the exponents one should go beyond mean-field, by including diffusion terms and noise in Eq. (|9|) 
and applying dynamic renormalization group techniques |T^, 
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